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Abstract 

It has recently been argued that non-trivial Brans-Dicke black hole solutions 
different from the usual Schwarzschild solution could exit. We attempt here 
to “censor” these non-trivial Brans-Dicke black hole solutions by examining 
their thermodynamic properties. Quantities like Hawking temperature and 
entropy of the black holes are computed. Analysis of the behaviors of these 
thermodynamic quantities appears to show that even in Brans-Dicke gravity, 
the usual Schwarzschild spacetime turns out to be the only physically relevant 

uncharged static black hole solution. 
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1. Introduction 


Brans-Dicke (BD) gravity [1] is perhaps the most well-known alternative theory of clas¬ 
sical gravity to Einstein’s general relativity. This theory can be regarded as an economic 
modification of general relativity which accomodates both Mach’s principle and Dirac’s large 
number hypothesis as new ingredients. Ever since it Erst appeared, it has remained as a vi¬ 
able theory of classical gravity in that it passed all the available observational/experimental 
tests provided a certain restriction on the generic parameter, “cn” of the theory is imposed 
[4], Shortly after the appearance of their first work [1], one of the authors, C. Brans pro¬ 
vided static, spherically-symmetric metric and scalar held solutions to the vacuum BD held 
equations [2], Since the gravitational collapse and the subsequent black hole formation is 
generally of great interest in classical gravity, in the present work we would like to address 
questions like ; under what circumstances Brans’ solutions can describe black hole space- 
times and if they actually do, they could really be non-trivial ones different from the general 
relativistic black holes ? It will then be discussed that although non-trivial BD black hole 
solutions different from the usual Schwarzschild solution appears to exit as suggested re¬ 
cently by Campanelli and Lousto (CL) [3], when “censored” by quantum aspects of black 
holes, namely their thermodynamics, it can be shown that they cannot really arise in nature. 
Brans [2] actually has provided exact static and isotropic solution to the vacuum BD held 
equations in four possible forms depending on the values of the arbitrary constants appear¬ 
ing in the solution. In the present work, we consider only the Brans ‘type I' solution since it 
is the only form that is permitted for all values of the “BD parameter”, uo (the other three 
forms are allowed only for negative values of uj, uj < —3/2). In fact in this work, we will 
exclusively assume that the parameter uo is positive since it has been prescribed so originally 
in the BD theory itself [1] (namely according to Brans and Dicke, the positive contribution 
of nearby matter to the spacetime-dependent Newton’s constant demands uj be positive) and 
it also has been constrained so by experiments [4] (for uj > 500, the theory is in reasonable 
accord with all available experiments thus far.) In fact, there is another crucial reason why 
uj has to be positive from held theory’s viewpoint. Namely, in order for the BD scalar held 
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$ to have “canonical (positive-definite)” kinetic energy, u needs to be positive. 

The vacuum BD gravity is described by the action (we shall work in the Misner-Thorne- 
Wheelcr sign convention) 

s = wJ ] (i) 

and the field equations 

R,,, - \g^R = - ^V a $V“*] 

+ ^ [V M '^ V Q V"$], (2) 

V Q V“<h = 0. 


The Brans type I solution obtained in isotropic coordinates is given by [2] 

4 / ^ 2(1 — Q) 

r 0 ' 


ds 2 = 


r — \ 2 (Q~X) 

r —^) cte 

r + r 0 / 


4 + 

$(r) = Q 


' r — tq ^ x 


+ r 0 


With g 2 + (l + |)y 2 -gy-l=0. 


( 3 ) 

( 4 ) 


Here, the arbitrary constants Q : x appearing in the solution are subject to the constraint in 
eq.(4) and they are related to Brans’ original notation [2] by Q — (1 + c)/A, y = c/A and 
to the notation of CL by Q = (1 — n), x — ~( m + n )- Next, $(r) denotes the “Brans-Dicke 
scalar field” and the quantity To, which is related to Brans’ notation [2] by r 0 = B (or by 
r 0 = ro to the notation of CL), is a “mass parameter” related to “scalar” and “tensor” mass 
by M s = — yr'o, M t = (2 Q — y)r 0 respectively. The quantities M s and M t are constants 
and related to the “Keplerian” (active gravitational) mass M measured by a test particle 
by M — M t + M s . Like the ADM mass in general relativity, the tensor mass M t is positive- 
definite and decreases monotonically by emitting gravitational radiation. The scalar mass 
M s and the Kepler mass M, however, share none of these properties [5]. Also note that for 
parameter values Q = 1, y = 0 (or m — n — 0 in the notation of CL), the above exact 
solution in eq.(3) reduces to the usual Schwarzschild soluton in Einstein gravity implying 
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that it indeed is a particular solution of BD field equations. 

Now, consider the coordinate transformation to the standard Schwarzschild coordinates 

r = r(l H-) 2 or f = -[(r — 2r 0 ) + (r 2 — 4-ror) 1 / 2 ]. 

f 2 

In terms of Schwarzschild coordinates, then, the Brans type I solution in eq.(3) now takes 
the form 

ds 2 - -(1 - —)^~^dt 2 + _ - _r/r 2 4-_ — _ dVt 2 

aS - r > at + ( 1 _4ro )Q ar + (1 _ 4roj(Q-l) fl ^ 2 ’ 

$(r) = (l-^)f. (5) 

r 

2. Non-trivial classical BD black hole solutions 

As mentioned earlier, CL recently examined this vacuum solution and claimed that under 
certain circumstances it could represent non-trivial BD black hole solutions different from 
the usual Schwarzschild solution. In what follows, we first give a brief review of their 
argument. Obviously in order for this metric solution to represent a black hole spacetime, 
the parameters ( Q , y) appearing in the solution should satisfy certain conditions. One 
straightforward way of obtaining such conditions is to consider under what circumstances 
an event horizon forms. To do so, CL studied outgoing null geodesics and looked for a 
condition under which the surface at r = 4ro (or in isotropic coordinates, at r = r o) could 
behave as an event horizon. It turned out that it happens provided 

2Q- X >1. 

A simple argument that can lead to this condition goes as follows ; first, since the metric 
solution of the vacuum BD held equations above has time translational (t —> t + St) isometry, 
a timelike Killing held C' = exits correspondingly. Now, as is well known, in general if 
is a hypersurface orthogonal Killing held which commutes with all others (if there are 
more than one), then the surfaces where = 0 are “Killing horizons”. Thus for the case 
at hand, in order to see if a Killing horizon such as an event horizon develops, we need to 
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find out under what condition can have zeros. It then is straightforward to see that 
£ /J £ /t = g 00 = (1 — 4r 0 /r)^ -x ) = 0 can have a zero at r = 4r 0 provided (Q — x) > 0- I n 
addition, certainly we do not want to have a curvature singularity at the surface r = 4ro 
which is a candidate for an event horizon. Thus from = r 2 (l — 4ro/r) (4_< 3) sin 2 9, we 
demand (Q — 1) > 1. These two conditions, when properly put together, then yields the 
above condition 2 Q — y > 1 which may be thought of as the condition for this metric solution 
in eq.(5) to represent possibly a black hole spacetime. 

Next, in order for this metric solution to represent truely a black hole spacetime, it should 
have a “regular” event horizon. As mentioned briefly above, then, we should further require 
that the curvature, for instance, have a non-singular behavior at the null surface r = 4r 0 . 
Thus in this time we consider the Kretschmann curvature invariant ; 

I = ( 6 ) 

= - ln ) P«-4) {( 2n ) 2 /iWi , + 4 (^) i 2 ( q , X ) + 6 i 3 ( q , x)} 

y*\J y* y* y* 

where 

h(Q, x) = 7 Q 4 + 16Q 3 + 14 Q 2 + 8 Q + x 4 + 2y 3 + 6y 2 - 16g 3 y 
+ 15Q V - 6 Q X 3 ~ 51 Q 2 X + 20 Q X 2 - 12 Q X + 3, 
h{Q, x) = -4Q 3 - 6 Q 2 + 11Q + y 3 - 3y 2 + 7Q 2 y - 5Qy 2 + QQx ~ 13, 
h(Q j x) — 2 Q 2 + x 2 ~ 2 Qx- 

Note hrst that this curvature invariant goes to zero as r —> oo as fast as / —» 0(r -6 ) and 
for the special case of interest, Q = 1, y = 0, i.e., for the Schwarzschild solution, it reduces 
to I = 48 ^6°' 1 = 48A S /2 (where M = 2r 0 is the ADM mass for Schwarzschild solution) as it 
should. Next, it is easy to see that the condition for non-singular behavior of the curvature 
invariant at r = 4r 0 amounts to the constraint 

Q >2. 

Finally, put them altogether, the condition for the metric solution in eq.(5) to represent a 
black hole spacetime with regular event horizon turns out to be 
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2 Q — x > 1 and Q > 2. 


( 7 ) 


Therefore it appears that for these values of the parameters appearing in the solution, 
the metric in eq.(5) could represent a non-trivial black hole spacetime different from the 
usual Schwarzschild solution in general relativity. At this point let us recall the well-known 
Hawking’s theorem on black holes in BD theory [11]. Long ago, Hawking put forward a 
theorem which states that stationary black holes in BD theory are identical to those in 
general relativity. To be more concrete, Hawking extended some of his theorems for general 
relativistic black holes to BD theory and showed that any object collapsing to a black hole 
in BD gravity must settle into final equilibrium state which is either Schwarzschild or Kerr 
spacetime. And in doing so, he assumed that the BD scalar field $ satisfies the weak energy 
condition and is constant outside the black hole. Now one may be puzzled. It is being 
claimed that non-trivial BD black hole solutions different from the Schwarzschild black hole 
might exit in apparent contradiction to the Hawking’s theorem. The possible existence 
of non-trivial BD black hole solutions discussed here, however, does not really contradict 
Hawking’s theorem since they violate the weak energy condition on the BD scalar field <f> 
that we now show below. 

The energy-momentum tensor for the BD scalar field is given by 

(*) = A[p(V„*V„* - + iv„V F 4>] (8) 

from which one can readily compute 

= Too = Q 2 - 1)(—) 2 (1 - — ) 2(Q - 1} -*. (9) 

o7r r r 

Thus for Q > 2, which is the condition satisfied by regular BD black hole solutions, 
< 0, namely the weak energy condition on the BD scalar field is violated. And 
this is the only means for the non-trivial BD black hole solutions to evade Hawking’s the¬ 
orem. Since they violate weak energy condition on BD scalar field, one may simply reject 
them as unphysical solutions. One may, however, still take them seriously as CL did by 
keeping the viewpoint that demanding the weak energy condition on the BD scalar field is 
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not absolutely compelling. In the present work, we choose to take the viewpoint of the latter 
and then proceed further. In the next section, we shall attempt to “censor” these non-trivial 
BD black hole solutions by examining their thermodynamic properties. Our philosophy here 
is that not all classically allowed black holes might be truely realistic and they should be 
further censored by quantum aspects of black holes, i.e., their thermodynamics. 

3. Thermodynamics of the non-trivial BD black holes 
We now turn to the investigation of the thermodynamics of non-trivial BD black hole solu¬ 
tions discussed thus far. The long-standing gap between thermodynamics and gravity has 
been essentially bridged by Hawking and later by many authors [6] who first have shown, via 
the study of quantum fields propagating on black hole background spacetimes, that black 
holes do evaporate as if they were black bodies. Namely the black hole thermodynamics is 
essentially associated with the quantum aspect of the black hole physics. Therefore one may 
expect that by examining thermodynamic properties of a black hole, one can, to some extent, 
explore its quantum aspect. Besides, since the practical study of black hole thermodynam¬ 
ics begins and ends with the evaluation of black hole’s temperature and entropy, we shall 
attempt to compute them. First, in order to obtain the Hawking temperature T H measured 
by an observer in the asymptotic region, one needs to compute the surface gravity k of the 
black hole and relate it to the temperature by Th = h/2tc [6]. Although the identification 
of the black hole temperature with Th = n/2i r was originally derived from studies of linear, 
free quantized fields propagating in given black hole geometries [6], it holds equally well for 
interacting fields of arbitrary spin and for general black hole spacetimes [9]. Thus our task 
reduces to the calculation of the surface gravity. In physical terms, the surface gravity k is 
the force that must be exerted to hold a unit test mass at the horizon and it is given in a 
simple formula as [10] 

k 2 = 4(V'V)(v„x„) (io) 

= -j-fVavrW 

where is a Killing held of the given stationary black hole which is normal to the horizon 
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and here the evaluation on the horizon is assumed. For our static BD black hole given in 
eq.(5), is just the timelike Killing field, and hence the surface gravity and the 

Hawking temperature are computed to be 

t « = T = T (q - xX^ttX 1 - (11) 

2n 2ir r z r 

Having established the expression for the Hawking temperature X# of the black hole, we 
next go on to find the “local temperature”. Generally, the local (or “Tolman redshifted”) 
temperature T(r) of an accelerated observer can be obtained by blueshifting the Hawking 
temperature Th of an observer in the asymptotic region from infinity to a finite point r. In 
other words, according to “Tolman relation” [7], the local temperature T(r) measured by 
a moving, accelerated detector is related to the Hawking temperature T H measured by a 
detector in the asymtotic region by 

T{r) = (-g 0 or 1/2 T H - (12) 

Thus the local temperature of our BD black hole is found to be 

T(V) = 2-(Q - xU^rUl - — ) ,Q - x/2 - ,, lr.4 n x (1 - ( 13 ) 

which behaves asymptotically as T(r —> oo) —> Th as it should. Next we turn to the 
computation of the black hole’s entropy. Generally speaking, the black hole entropy can be 
evaluated in three ways ; firstly, following Bekenstein-Hawking proposal [8], one can argue a 
priori that the entropy of a black hole must be proportional to the surface area of its event 
horizon (i.e., S = \A). Alternatively, knowing the Hawking temperature Th and chemical 
potentials (i.e., coulomb potential at the event horizon for the conserved U(l) charge Q 
and angular velocity of the horizon Qh f° r the conserved angular momentum J generally), 
one may integrate the 1st law of black hole thermodynamics TndS = dM — ndQ — TtndJ 
to obtain the entropy [9]. Thirdly, according to Gibbons and Hawking [9], thermodynamic 
functions including the black hole entropy can be computed directly from the saddle point 
approximation to the gravitational partition function (namely the generating functional 



analytically continued to the Euclidean spacetime). In the present case where we consider the 
static, isotropic black holes in BD gravity theory, the expression for the Hawking temperature 
given in eq.(ll) renders it awkward to employ the second method which uses the 1st law 
of black hole thermodynamics. Nor can we naively adopt the Bekenstein-Hawking relation 
S = jA to obtain black hole’s entropy. In fact, the relation S = jA has been established 
essentially in the context of the Einstein gravity via the 1st law of black hole thermodynamics 
which has been derived from the expression for the Bondi mass of a stationary black hole 
[10]. In the context of BD gravity, clue to the addition of the gravitational scalar degree 
of freedom (i.e., the BD scalar field $), the expression for the mass of a hole is subject 
to a modification which, in turn, would result in another modification in the 1st law of 
thermodynamics violating the exact relation, S = \A. Therefore, here, as a reliable way of 
evaluating the entropy, we should resort to the method suggested by Gibbons and Hawking 
[9] since it is based on the fundamental definition of entropy. To illustrate the procedure 
briefly, we begin with the vacuum Euclidean BD gravity action and field equations 

I[g,$\ = - J- [| y d 4 x - ug^ + 2 d 3 xVh$(K - if 0 )] • (14) 

The boundary term in Euclidean BD action I[g, <h] above has been determined as follows ; 
following Gibbons and Hawking [9], we start with the form f dY cPx\/7i,B where B = — + 

C. Now for asymptotically flat spacetimes where the boundary dY of the 4-dim. manifold 
Y can be taken to be the product of the (analytically-continued) time axis with a 2-sphere 
of large radius, i.e., dY = S 1 x S 2 , it is natural to choose C so that B vanishes for the 
flat spacetime metric rj /w . Thus we have B = — K 0 ) where K and K 0 are traces 

of the second fundamental form of dY in the metric g^ u and rj^ w respectively (one may 
choose to take B = — gk(<h/l — $ 0 -^ 0 ) where $0 = constant = 1. The choice of the latter, 
however, leads to no substantial change in the conclusion that we shall draw from the 
former.). Now the gravitational partition function and the entropy are given in the saddle 
point approximation by 

z = fM* 
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(15) 


S = In Z + f3M ~ -I[g c , 4> c ] + — 

where the superscript “c” denotes the saddle point of the action. Since the saddle point of 
the action is nothing but the solution to the vacuum BD field equations in eq.(2), first we 
have 

J y d 4 x Jg (<M - = 0 

at ( g c , <f> c ) and thus the non-vanishing contribution to I[g c , < f >c ] comes only from the boundary 
term on dY. Here the computation of the term f 9Y d 3 x\fh<&KQ is straightforward since we 
know that for flat spacetime, Kq = 2 jr. Next from the definition of the second fundamental 
form K [10] and from the spherical symmetry of the BD black hole solution given in eq.(5), 
we have f 9Y d 3 x\fh<&K = <I>^- f 9Y d 3 x\/h, where f 9Y d 3 x\/h, is the derivative of the area 
f 9Y d 3 x\fh, of dY' as each point of dY is moved an equal distance along the outward unit 
normal n. Then the result of the actual calculation is 

I[g c , <f> c ] = 27rr 0 fiT 1 (3<2 + y - 2) + O^r _1 ). 

Therefore the black hole entropy is found to be, 

S=[2-(Q + 3x)]^ (16) 

d-H 

_ . 2 - (<3 + 3x) , 4r 0'-(O- > /2-l) 

“ 1 (Q-x) 1 1 r ' ' r ~ ir °' 

where we have used M = M t + M s = 2 r 0 (Q — y). Note that the entropy S is inversely 
proportional to the temperature Th and it does not satisfy the usual Bekenstein-Hawking 
relation [8], S = - { A = 7rr 2 (l — ^ m ) _ ^ _1 )| r=4ro as speculated earlier. Although we have 
obtained the Hawking temperature T Hl the local temperature T(r) and the entropy S of 
our BD black hole in eq.(5), no dehnite statement concerning their behaviors can be made 
unless the precise estimation of the allowed values of the combination of parameters such 
as ( Q — y/2) and ( Q — y) are done. Therefore, next we carry out the estimation of allowed 
values of these parameters. In determining the allowed values of them, two bottomline 
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conditions are ; (i) the constraint equation in (4) that parameters Q and y should satisfy 
must always be imposed and (ii) the condition for the possible formation of non-trivial BD 
black holes given in eq.(7) must hold. First, as for the quantity ( Q — y/2), the condition (ii) 
gives (Q — y/ 2 ) > 1/2 and Q > 2 whereas the condition (i) leads to — 1 < (Q — y/ 2 ) < 1 
which, when put together, yields 

\ < (Q - f) < 1 and Q > 2 . (17) 

Next, it is now the turn of the quantity (Q~x) an d the condition (i) leads to < ( Q — 

x) < 3 while the condition (ii) or equivalently the condition obtained in eq.(17) gives in 

addition to g > 2 , {(g-y) < |^±|}n{(g-y) < |( 1 -y^J^) or {Q~x) > + 

which, when put together, yield 


1 2a; + 4 1 

-< (g — y) < -(i — 

2u + 3 J 2 v 


2u + 3 ' 


and g > 2 


( 18 ) 


or 

5 (1 + K iQ ~ x) ~ /HI alld °- 2 ' 

Further, from the relationship between a = (g — y/2) and (3 = (Q — x)i that is, 
ck = (2ca + 4) _1 [(2cu + 3)/? ± (2cu + 4) — (2cu + 3)/3 2 ], it is straightforward to see that, 
for parameter values obtained in eqs.(17) and (18), the overlap 1/2 < a < 1 and 
1(1 + y 2 ^+ 3 ) < P < can happen as is manifest, for instance, for a = l/\/2, 

f3 = ^ 75(1 + ^ 2 ^+ 3 ) w h erea s the other overlap can never happen since we assume u> > 0. 
In addition, consider as cases of special interest, g = 1, y = 0 and g = 1, y = -/h. The 
first case is the Schwarzschild solution as mentioned earlier and the second case corresponds 
to a BD black hole with a singular horizon but satisfying the weak energy condition on the 
BD scalar held. Finally, consider three cases (including two cases of special interest just 
discussed) 

case (I) 1/2 < (Q - y/ 2 ) < 1 , \(l + y/JQ C (Q - y) < ylSl and Q - 2 

case (II) g = 1 and y = ^7 
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case (III) Q = 1 and x = 0 

The behaviors of thermodynamic functions T H , T(r) and S in eqs.(ll), (13) and (16) re¬ 
spectively are given for each of above three cases as 
case (I) Th = oo, T(r) = oo, S = 0 
case (II) T H = oo, T{r) = oo, S = 0 

case (III) T H = j±r o , T(r) = S = 16nr 2 0 

This behavior of thermodynamic functions reveals that BD black hole solutions with pa¬ 
rameter values other than Q — 1 and x = 0 appears to be physically irrelevant and hence 
should be rejected as unphysical solutions. Obviously this is based on the observation that 
for cases (I) and (II), we have S = 0 and T H — oo at any stage of the Hawking evaporation 
which is against our conventional wisdom based on the usual statistical mechanics. One 
might want to propose different interpretation of these behaviors of thermodynamic quan¬ 
tities. Namely, since an infinite value of temperature might signal the breakdown of the 
semiclassical treatment involved in typical black hole thermodynamics, it may well be that 
the non-trivial BD black holes (?) could already be quantum entities which do not admit 
semiclassical treatment. It appears that the breakdown of the semiclassical approximations 
can be attributed either to the failure of weak energy condition or to the singular behavior 
of the horizon. Indeed the case (I) corresponds to regular black holes violating weak energy 
condition on whereas the case (II) represents a singular black hole satisfying weak energy 
condition. Generally, it seems that the non-singular behavior (which requires Q > 2) and 
the weak energy condition on $ (which demands —1 < Q < 1) cannot be simultaneously 
accomodated by any non-trivial BD black hole. However, since we can always retain non¬ 
singular behavior at the expense of weak energy condition (simply by demanding Q > 2), 
the breakdown can be attributed solely to the latter. Thus to summarize, non-trivial BD 
black hole solutions violating the weak energy condition on $, that we allowed classically, 
turned out to be entities which demand more rigorous quantum treatment. Now we seem 
to have two options ; first, semiclassical analysis in the black hole thermodynamics indeed 
works and hence the non-trivial BD black holes should be rejected as they fail this quantum 
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censorship. Second, they are really quantum entities to which the conventional black hole 
thermodynamics cannot be naively applied in the first place. In the absence of a consistent 
theory of quantum gravity, the second option is well beyond our scope and in this option we 
do not even know whether or not these non-trivial BD solutions can be identified with black 
holes at all. Thus in the present work, we choose to take a more conservative viewpoint, 
namely the first option. 

4. Discussions 

To conclude, in the present work we examined the thermodynamic properties of all classi¬ 
cally allowed static, spherically-symmetric black hole solutions in vacuum BD theory. And 
the lesson we learned from our analysis is as follows ; not all classically allowed black hole 
solutions may be physically realistic. They should be further censored by quantum nature 
of black holes, namely their Hawking evaporation mechanism. And when applied to a par¬ 
ticular case discussed in the present work, it turned out that the non-trivial solutions that 
violate the weak energy condition on the BD scalar field and thus have been discounted in 
Hawking’s theorem fail to survive the quantum censorship and hence would not really arise 
in nature. After all, stationary black holes in BD theory appear to be identical to those in 
general relativity once they settle down. 
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